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We prove the Razumov-Stroganov conjecture relating ground state of the 0(1) loop model 
and counting of Fully Packed Loops in the case of certain types of link patterns. The main 
focus is on link patterns with three series of nested arches, for which we use as key ingredient 
of the proof a generalization of the Mac Mahon formula for the number of plane partitions 
which includes three series of parameters. 
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1. Introduction 

The Razumov-Stroganov (RS) conjecture [1] relates the components of the ground 
state of the 0(1) loop model, which are indexed by link patterns (pairing of points on a 
circle) to the numbers of Fully Packed Loop configurations (FPL) on a square grid with 
a connectivity of external vertices given by the link pattern. Despite considerable activity 
around this conjecture 0||[||||6|,[F],|8| , it has not been proved yet. It is the author's belief, 
however, that the work || was a significant step in this direction: in it, an inhomogeneous 
loop model was introduced in order to make the ground state a polynomial of the inho- 
mogeneities (spectral parameters). This way, a corollary of the RS conjecture (which was 
already formulated in [|IU]]), namely that the properly normalized sum of all components 
of the loop model ground state equals the total number of FPL, also known as the number 
of Alternating Sign Matrices, was proved. 

The present work tries to demonstrate in a very simple setting how the methods of 
could help to prove the RS conjecture by considering a special subset of possible link 
patterns, namely those with few "little arches" (arches connecting neighbors). Consider 
the link patterns of Fig. 1. They are made of three sets of a, b, c nested arches. Here a, 
6, c are arbitrary integers such that a + b + c = n where the size of the system is In. The 
model depends on 2n complex numbers a±, . . . , a&+ c , Pi, ■ ■ ■ , Pa+c, 7i> • • • > 7a+b which are, 
up to multiplication by a power of q (as will be explained below), the spectral parameters 
of the model. Note there are several good reasons to restrict oneself to such link patterns, 
a particularly obvious one being that we know the corresponding number of FPL: it was 
computed in - and happens to be equal to the number of Plane Partitions in a hexagon 
of shape a x b x c! 




c+1 

Fig. 1: Link pattern with 3 sets of nested arches. 
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The reader is referred to 0] for details concerning the 0(1) loop model. We briefly 
recall its definition here, if only to fix notations. Let n be an integer, and consider a 
complex vector space equipped with a basis indexed by non- crossing link patterns: the 
latter are by definition pairings of 2n points on a circle, in such a way that the pairings 
can be represented by non-crossing edges inside the circle. On this space we define the 
action of a linear operator, the so-called transfer matrix T n (t\z\, . . . , Z2n) 5 which depends 
on complex parameters t,Zi, . . . , Z2m by the following graphical description: 

TMzu .... Z2n) = ft (f^Qa + i^a) (1.D 

where q = e 2l7r / 3 , and the symbolic product over i means that the plaquette of index 
i should be inserted at vertex i. The result is that, starting from a given link pattern, 
the action of T n produces a new link pattern by adding a circular strip of plaquettes and 
removing any closed loops thus created. The coefficients in Eq. ( |1.1| ), in the range of 
parameters where they are real, are simply probabilities of inserting the corresponding 
plaquettes, parametrized in a convenient way in terms of the Zi/t. As a consequence of 
the Yang-Baxter equation, [T n (t), T n (t')] = with all other parameters Zj fixed. 

Since T n is a stochastic matrix, it has the obvious left eigenvector (1,...,1) with 
eigenvalue 1. Therefore, it also has a corresponding right eigenvector, which is unique for 
generic values of the zf. 

T n (t\Zl, . . . , Z2n)^n(z 1 , . . . , Z 2n ) = • • • , Z 2n ) (1.2) 

One can normalize \I/ n in such a way that its components in the basis of link patterns 
7r are coprime polynomials of the Zi. One still has an arbitrary numerical constant in 
the normalization of \E' n . Consider now the homogeneous limit when all Zi equal 1. If 
one chooses this constant so that the smallest entry is 1, then the \E' n)7r (l, . . . , 1) are the 
subject of various conjectures, including the remarkable Razumov-Stroganov conjecture 
already mentioned above, that identifies them with a certain FPL enumeration problem. 
In what follows we shall choose another numerical normalization which is more convenient 
for intermediate calculations. 

In Sect. 2, we derive the main formula for the entries of the ground state of the 0(1) 
loop model corresponding to link patterns with three sets of nested arches (Fig. 1). In 
Sect. 3 we establish the connection with plane partitions (or dimers). Sect. 4 discusses 
the (partial or total) homogeneous limit. Sect. 5 briefly describes the extension to more 
general link patterns for which the corresponding enumeration of FPL is known. Sect. 6 
concludes. 
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2. Recurrence relations and their solution 

In ||, a certain number of relations were shown to be satisfied by \I/ n , the ground state 
eigenvector of the 0(1) loop model of size 2n. We shall need the following three facts: 

Theorem 4 of H . The components of ty n are homogeneous polynomials of total degree 
n(n — 1), and of partial degree at most n — 1 in each variable Z{. 

Theorem 1 of |J . The entries \I/ ni7r of the groundstate eigenvector satisfy: 

*n,7r(^l, • • -,Z 2 n) = I ] [ ] [ i'Fi ~ ' •-./ ) (~1 =2» ) (2.1) 

\ s€E T i,jEs / 



where $ ni7r is a polynomial which is symmetric in the set of variables {z{, i6s} for each 
s G E n , and E n is the partition of {1, ... , 2n} into maximal sequences of consecutive points 
not connected to each other by arches of n. 

Theorem 3 OF ||. If two neighboring parameters z% and Zi+i are such that Zi+i = q 2 Zi, 
then either of the two following situations occur for the components ^ n ^: 

(i) the pattern it has no arch joining i to i + 1, in which case according to Theorem 1, 

*»,7r(zi, Zi+i = q 2 Zi, . . . , z 2n ) = ; (2.2) 

(ii) the pattern tt has a little arch joining i to i + 1, in which case 



^n,7r(zi, . . • , Zi, Zi+l — q 2 Zi, . . . , Z2n) 
( ^ \ 



n ) (2-3) 



where it' is the link pattern 7r with the little arch i, i + 1 removed. 

In what follows we shall concentrate on components corresponding to link patterns 
with three sets of nested arches of size a, 6, c, which we shall denote by ^ a ,b,c- The spectral 
parameters are relabelled as Z\ = ctj, qfy, q 2 ji according to the pattern of Fig. 1. Thanks 
to Theorem 1, we can write 

y a , b ,c= n (9« l -9" i « J ) n (qfc-q' 1 ^) n (^-^so^c^) 

l<i<j<b+c l<i<j<a+c l<i<j<a+b 
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where the arguments {zi} = {ctj, g/3j, q 2, ji} have been suppressed for brevity. According 
to Theorems 1 and 4, $ a ,6,c is a polynomial of total degree ab + bc + ca, and a symmetric 
polynomial of the {ttj} of degree at most a in each, of the {(3i] of degree at most b in each, 
and of the {7^} of degree at most c in each. 

We now rewrite Eq. ( |2.3| ) in the case when z, is the last parameter a and Zi+i is the 
first parameter /3, in terms of $ a ,6,c- Since the latter is a symmetric function it is actually 
irrelevant which a and which (3 are singled out, and the result is: 

a+b 

$o,6,c|/?j=a« = Y[(®i - 7k) $a,6 >C -l (2.5) 
fc=l 

where the parameters of <& a ,6,c-i are the same as those of $ a ,&,c, except aij and are 
removed. Since <& a ,b,c is of degree a in a„ the equations ( |2.5|) with j = 1, . . . , a + c and 
fixed i determine entirely $ a ,6,c as long as c > 1. They form a very simple recurrence 
relation which is supplemented by the initial condition $ a ,6,o : this corresponds to the 
so-called "base link pattern" (Fig. 2), which is entirely factorized by Theorem 1: 

b a 

*40 = nil(«*-ft) ( 2 - 6 ) 

i=lj=l 




a+b 



Fig. 2: Base link pattern. 



We now claim that the following Ansatz solves the recurrence relation: 

ra+c/ n \ TT T-ra+6 . 



E 



IC{l,...,b+c} 

#I=c 



nieillj0j(a*-«i) 



(2.7) 



The case c = forces 7 = and we recover immediately Eq. (|2.6| ). Next, notice that 
the expression of Eq. ( |2.7|) is symmetric in all three sets of variables: it is obvious for the 



{(3i) and the {7^}; it is also clear for the {o^} since the summation over all possible subsets 
of cardinality c is invariant by permutation of {1, ...,&+ c}. We can therefore choose one 
CKj and one (3j and set them equal, say (3 a+c = a^+c- This forces b + c G I in Eq. ( |2.7|) . 
Define J' = 7 — {6 + c}. The summation over / can then be replaced with the summation 
over /' C { 1 , . . . , 6 + c — 1}, and it is easy to check that cancellations in numerator and 
denominator reproduce Eq. ( |2.7| ) with c — > c — 1 and the parameters otb+c, /3 a + c removed. 
This concludes the recurrence. 

Note immediately the symmetry between the sets of variables {/3i} and {7^} in 
Eq. ( |2.7| ): indeed, replacing / with its complement I exchanges their roles (as well as 
b and c). However, the {cti} seem to play a different role. We shall now produce an 
equivalent expression which restores the symmetry a <-> /3, at the expense of breaking the 
symmetry [3 <-» 7: 

b+c a+c 



^4nnK-ft-)/-/|-| n 

nLinS(^-Tfe) 



il=i njsc* - nLi n-i c (^ - 



The c contour integrals should be defined in such a way as to encircle (counterclockwise) 
all the poles ol\ (but none of the Pi). One goes back to Eq. ( |2.7|) by applying the Cauchy 
formula. Each z% must be evaluated at a certain ai i with 1 < Ii < b + c; furthermore, the 
factors Yl( z i ~ z j) force the Ii to be distinct, and we reproduce after various cancellations 
the summation over I = {ii, . . . , h+c} of Eq. (|2.7|) . 

The formula ( |2.8| ) is of the form of a matrix integral: the contour integral makes it 
essentially similar to the unitary matrix integral. This analogy will be pursued below. 
For now, we use a standard trick in random matrix theory, which is to introduce the 
Vandermonde determinant A(zj) = n^j^i ~ z j) = det(z^ _1 )i<i< c , and then to note that 
det (Pi (zj)) = A(zi) detP where the Pi, 1 < i < c, are arbitrary polynomials of degree less 
than c and P is the c x c matrix of coefficients of the Pj. In Eq. ( |2.8|) we have a squared 
Vandermonde determinant, so we can introduce another similar set of polynomials Qi. 
Moving the determinants out of the integrals, we find: 

nSiE2(°*-& 



a ' b ' c det P det Q ^ 



<f P £ (z)Q m (z) 
J 2m 



dz (~\ Y[lt\(z ~ Ik) 

ra+c / n \ TT^+c / 



In what follows, we shall be naturally led to a choice of polynomials P and Q. 



l<£,m<c 

(2.9) 
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3. Connection with plane partitions 

We now introduce a model of weighted Plane Partitions - in more physical terms, it 
is a model of dimers on the hexagonal lattice, but we shall prefer the language of Plane 
Partitions in what follows. Configurations are defined as tilings with lozenges of a hexagon 
of size a x b x c. Lozenges come in three orientations since they are made of two adjacent 
equilateral triangles of a regular triangular lattice. The model comes with three series 
of parameters ai, [3j, 7^ living on the lines of the underlying medial Kagome lattice, 
see Fig. 3 (i). To each lozenge of the plane partition (or equivalently to each dimer) is 
associated a local Boltzmann weight ctj — (3j, 7^ — (3j, ai — 7^ given by the difference of 
the parameters of the lines crossing at its center. Note that there are exactly ab, be, ca 
lozenges of each orientation. 




Fig. 3: (i) Plane partition and its parameters, a = 2, b = 4, c = 3. (ii) 
Associated non-intersecting paths. 

We wish to compute the partition function Z a ^ jC , i.e. the sum over all configurations 
of the product of local Boltzmann weights. In order to do so, it is convenient to use yet 
another representation in terms of non-intersecting paths (NIPs). To each plane partition 
one can associate c paths originating from one of the sides of length c and ending at the 
other, see Fig. 3 (ii), which simply follow two types of tiles out of the three. 

One can replace the local Boltzmann weights of the tiles with a local probability for 
the path to go left or right, by factoring out all the possible weights of the third type of 
tile: 

Z a ,b,c= j I (Cti-/3j) F a;b . c (3.1) 

l<i<b+c,l<j<a+c 
i+j>c, i+j<a+b+c 
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and then introducing an inverse weight when the tiles of the third type are absent (i.e. 
where the paths are): 

{— - — edge at the crossing of (ctj,7fc) 
T-% ' 3 - 2 > 
-f- edge at the crossing of (/L, 7*.) 
OLi - Pj ~ 

where the two orientations of the edges (or of the underlying dimers, or lozenges) determine 
which weight to use, and the third coordinate is given by i + j = k + c. These NIPs move 
exactly o steps in one direction and b steps in the other. 

NIPs are free fermions, and therefore their propagator F a ^. c is a determinant of one- 
particle propagators: 

F a ,b;c = det [<^moii (a, 6, c)] 1 < £m < c (3.3) 

where I — > m means the probability (for a single path) to go from position I on one side of 
length c to position m on the other. This is also known as the Lindstrom-Gessel-Viennot 



formula |T2]. The one-particle propagator is nothing but the case c = 1, with appropriately 
shifted a, 6, a,, j3y. 

[£ ^ m on (a, b, c)] 

Fa-\-£— m,b-\-m— £,1 

(a>£, . . . , a b+m ; /3 c +i-i, ■ ■ ■ , Pa+c+i-m", 7i> • • • ; 7a+b)) 

(3.4) 

We are finally led to a simple problem of computing the weighted enumeration of a 
single path. The following formula holds: 

Znt U i= i {z - on) lii=i (* - Pj) 
(3.5) 

where once again the contour integral encircles clockwise the on but not the f3{. This can 
be proved by noting that -F a ,6;i satisfies the following simple recurrence formula: 

a b+1 - 7 a+b 7 a+b - (3 a+ i 

*a,b;l — a — ^a-l,6;l H 7, # a,b-l\\ ^.DJ 

«b+l — Pa CKf, — p a+ i 

and by using z - la+b = - &+0 + £5^fer(* " 

Putting together Eqs. ( |3. 1| ) (|5~5|) , one obtains 



Z a ,b,c= I I (oi-fij) det 



l<i<b+c,l<j<a+c 
i+j>c, i+j<a+b+c+l 



Ukt b i(z - Kk) 



l<£,m<c 

(3.7) 
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To connect with Eq. (|2.9|), set 

P e (z) = {z-ot 1 )---{z- ai-i) (z-Pi)~-[z- Pc-i) 
Qm(z) = (z - a b+m+1 ) ■■■(z- a b+c ) (z - (3 a+c+2 - m ) • • • (z - (3 b+c ) (3.8) 

By factor exhaustion one finds immediately that detP = rL+j<c( a i — Pj) an< ^ detQ = 
n i+j > a+b+c+2 (®i ~ Pj)- Plugging this into Eq. Q reproduces exactly Eq. ( gjg) . 

We conclude that Z a ^ b)C = <& a ,b,c- We have thus obtained a direct, exact relation 
between the components of the inhomogeneous 0(1) loop model corresponding to three 
sets of nested arches (a, 6, c) and the partition function of weighted plane partitions on a 
hexagon a x b x c. 



4. Relation to unitary matrix integrals and homogeneous limit 

In order to prepare the ground for the homogeneous limit, one can now choose some 
of the variables to be equal: on = a, Pj = f3. An important property of the matrix 
integral-like expression Q2.8| ) is that it is preserved by homographic transformations. We 
define w = |5§ to send (a, 0) to (0, oo) and obtain 



7fe-« % 



i<*<i< c 



C /" /dwi dw c -,-r ... _ u Ut=iWH\(wi 



7fe — Q! \ 

where C = (a — /3) a6 n&=i(7fc ~ PY- The are integrated on contours surrounding 0, for 
example |it^| = 1. We recognize in the second line of Eq. (|4.1|) the usual form of a matrix 
integral over the unitary group U(c) once angular variables are integrated out and only 
the eigenvalues Wi are left. We thus obtain 

$ abc = C' [ dfi det(l + T® 0)(detO _1 ) b (4.2) 

Ju(c) 

where dfi is the Haar measure on U(c), V is the (a + b) x (a + b) diagonal matrix with 
eigenvalues ] i k Z y ^ k ; an d C" = (a — PY b Y[k= b i(lk ~ a Y- O ne can now use the identity 
det(l + T ® O) = J2x s x(^)s\t(T), where A is a partition or Young diagram, s\ is the 
corresponding GL character (with the convention that it is zero if the Young diagram 
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has more rows than the size of the matrix), and A T is the transposed Young diagram. 
Orthogonality of characters of U(c) finally allows to perform the integration over Q and 
results in the simple formula 

$ a , b , c = C's Yb>c (T) (4.3) 

where Y"b jC is the rectangular Young diagram with b rows and c columns. Of course, this 
is a standard result once reinterpreted in terms of non-intersecting paths (sy bc (r), as a 
Schur function of the eigenvalues of T, can be defined as a sum over semi-standard Young 
tableaux, which are themselves in bijection with the NIPs). 

We finally consider the homogeneous situation where all Zi are equal to 1, that is 
Oii = a = 1, (3 j ; = (3 = q 2 , 7^ = 7 = q. In this case V = — q 1. We find that & a ,b,c becomes 
^(ab+bc+ca)/2 ^[ mes dimension of the GL(a + b) representation with rectangular Young 
diagram b x c. The latter is one of the many formulae for the number of plane partitions. 

The Razumov-Stroganov conjecture [1]] claims that in the homogeneous O(l) loop 
model, ty a ,b,c/^n,min (where ^ n ^ m in is the smallest component of \l/ n ) must be equal to the 
number of Fully Packed Loop configurations (FPL) with the corresponding connectivity 
(a, 6, c). With our normalization conventions, ^n,min = 3 n ( n_1 - ) / 2 and all powers of 3 
cancel out, so that ^ a ,b,c/^n,min is simply the number of plane partitions in the hexagon 
a x b x c. But according to [ I T|, the number of FPLs with connectivity (a, 6, c) is the very 
same number. This proves the RS conjecture for the case of these link patterns. 



5. Generalization to four little arches 

Let us first reobtain the result of the previous section in a more synthetic way. We 
use the fact, proved in appendix A, that when one switches the two spectral parameters 
of neighboring parallel lines, the partition function of plane partitions with an arbitrary 
geometry is unchanged. This implies that the partition function Z a ^^ c for plane partitions 
introduced in Sect. 3 is a symmetric function of the spectral parameters {o^}, {Pj}, {7fc}- 
At this stage, one can skip the entire reinterpretation in terms of free fermions and prove 
directly that it satisfies the same recurrence relations as $> a ,b,c- Indeed, setting say a.\ = 71 
forbids the lozenge parallel to sides a and c in the corner, see Fig. 4, and thus creates two 
rows of "frozen" lozenges which lead us back to the case a x b x (c — 1). This provides 
a nice graphical interpretation of the recurrence relations, very much in the spirit of the 
recurrence relations of Korepin for the six-vertex model with Domain Wall Boundary 



Conditions [13 
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a 



Fig. 4: Plane partition in which the a.\ and 71 rows are frozen. 




Z c+I 



Fig. 5: Link pattern with four little arches (a, b\e\c, d). 
Consider now the most general link pattern which possesses four "little arches" (arches 



connecting neighbors), as described by Fig. 5. It was shown in [[14[ that their enumeration 
is equivalent to that of certain lozenge tilings of a region of the plane with identifications, 



see Fig. 6. We refer the reader to [|15||14| for details. In particular note that there are 
exactly d "dents" in the two identified sides of length c + d. Inspired by the case of 
three little arches, it is natural to introduce spectral parameters into the lozenge tilings 
as described on Fig. 6. The weight of a lozenge is equal to qu — q~ 1 v where u and v are 
the spectral parameters crossing at the center of the lozenge in such a way that the line 
of v forms an angle of +n/3 with that of u (contrary to the case of three little arches, we 
cannot get rid of the factors of q by a redefinition of the spectral parameters) . One checks 
that this produces a partition function which has degree at most: c + d + e in each Xi, 
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Fig. 6: Lozenge tiling corresponding to the link pattern (3, 4|2|2, 1) and its 
parameterization. 

a + d in each t/i, a + b + e in each Zi, b + c in each £j, as should be. As a consequence of 
Appendix A, it is symmetric in each set of variables. 

It is now easy to check that these partition functions satisfy all the required recurrence 
relations. Among the various possibilities, two are depicted on Fig. 7. If one sets t\ = q 2 xi, 
the first rows of the sides of length a + d and c+d+e are frozen, and once these are removed 
one obtains the tiling with a — > a — 1. Similarly, if t\ = q~ 2 z±, the first rows of the sides 
of lengths b + c and a + d are frozen, and one dent becomes locked in first position. Once 
these rows are removed one recovers the tiling with d — ■> d — 1 (with, in particular, one less 
dent). 

Since the partition function is of degree b + c in t\ and is known at a + b + c + d 
values of ti, it is entirely fixed. The recurrence allows us to reach either d = or a = 0, 
at which point, up to an additional frozen region one recovers the tiling of a hexagon, a 
case which has already been treated. Now the component of \E' n corresponding to the link 
pattern (a, b\e\c, d), once rid of its factors Yii<j ( a x i ~Q~ lx j) Ili<j(iyi~ ( l~ 1 yj) Yl^jiQ^i — 
1~ lz j) Yli<j( a ^i ~ l^tj)) satisfies the very same recurrence relations (of course at each 
step one must check that the prefactors in the recurrence match); therefore it is equal to 
the partition function of lozenge tilings described above. In particular, when all spectral 
parameters are equal to 1, the homogeneous component is equal to the number of such 
tilings, as predicted by the RS conjecture. 
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Fig. 7: Lozenge tilings with (ti,Xi) and (ti,Zi) rows frozen. 
6. Conclusion 

In this article, we have proved in detail how the component ^ a ,b,c of the ground state 
of the 0(1) loop model corresponding to the link pattern with three series of nested arches 
(a, 6, c), once properly normalized, is equal to the number of plane partitions in a hexagon 
of size a x b x c. This is a highly non-trivial check of the Razumov-Stroganov conjecture 
since these link patterns form an infinite series. We have briefly described how the proof 
can be extended to more general link patterns, (a, b\e\c, d). 

In fact, we have found more than this: just as in H the sum of components was 
actually computed for arbitrary spectral parameters, here we have found that ^ a ,b,c with 
spectral parameters is the partition function of plane partitions with some local Boltzmann 
weights. This might seem unsurprising since one can hope that there is a unique "natural" 
way to introduce spectral parameters into the model; however, if ones tries to connect to 
FPL configurations, in which the RS conjecture is formulated, one is faced with a subtle 
problem: how to map back the plane partitions onto the square lattice of the FPL in such 
a way as to make sense of the spectral parameter dependence? Answering this is probably 
related to proving the full RS conjecture. We hope to come back to this point in the future. 

More obvious extensions of this work should be mentioned. Firstly, more general FPL 
configurations have been ]16[] or could in principle be enumerated: they have link patterns 
that can be obtained from those with four little arches by simple local modifications. 
Clearly, the formulae of || allow to obtain mechanically the corresponding ground state 
element. It would be a further check of the RS conjecture to take the homogeneous limit 
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in these formulae and recover the counting of FPL, though it is unlikely to produce any 
surprise. 

Secondly, An extension to arbitrary q of the polynomials was proposed in |17| and 



reformulated as a solution of the qKZ equation in |L8| . The present work being entirely 



based on recurrence relations which are still satisfied by solutions of the qKZ equation, it 
is clear that it can be extended to generic q. This might shed some light on a possible 
generalization of the RS conjecture to arbitrary q. 
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Appendix A. Proof of symmetry under interchange of spectral parameters 

The purpose of this appendix is to provide an explicit proof of the symmetry of the 
partition function of weighted plane partitions when one switches the spectral parameters 
of two neighboring parallel lines, say 71 and 72. The weighting is the same as described in 
Sec. 3 (i.e. for simplicity we absorb powers of q in the definition of the spectral parameters). 
The plane partitions are supposed to fill a domain with arbitrary shape; in fact, we shall 
only need to consider the region in which the two spectral parameters act and ignore the 
rest of the tiling. 

Generally, this region has the shape of Fig. 8. Note that the positions of the triangular 
holes are in principle not fixed: to a given shape of the entire domain to be filled may 
correspond several possibilities of locations for these holes. However here we shall keep 
them fixed and show that each term in the summation over them is a symmetric function 
of 71 and 72. The latter is equivalent to showing that a certain family of row-to-row 
transfer matrices is commutative; however here, to avoid any formalism we shall prove this 
by explicit computation. 




Fig. 8: Region of a tiling where the two spectral parameters 71, 72 appear. 
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The locations of the triangular holes are strongly constrained in order to allow the 
possibility of a lozenge tiling. As a consequence we see that we can divide our region into 
pieces separated by holes, so that these pieces have no influence on each other, and the 
summation over plane partitions factorizes accordingly. There are two types of pieces: 

a. Hexagons of size fe x 1 x 1, fe > 0. They are those that are delimited either (i) on at 
least one side by a pair of adjacent holes (by parity on the other side the boundary 
must have the same shape) or (ii) by two individual holes, with the parity being such 
that there are forced edges going inwards on both sides. The weighted enumeration 
of plane partitions inside a hexagon of size fcxlxl can be computed, cf Sec. 3: 

n?=iV» - 7i)(a - 72) - nSiVi - 72XA - 71) (A1) 

7i-72 

where the on and /?j are the spectral parameters in the other 2 directions. Eq. (|A.1| ) 
is explicitly symmetric in 71, 72. 

b. Parallelograms 2 x k, k > 1: they are those that are delimited by two individual holes 
on opposite boundaries, with the parity being such that there are forced edges going 
outwards. A parallelogram can only be filled by elementary lozenges in a unique way, 
hence the weight 

k 

n(a<-7i)(«<-72) (A.2) 

i=i 

where the oil are the spectral parameters going parallel to the non-horizontal sides of 
the parallelogram. 

We finally conclude that the whole partition function, being a sum (over the locations 
of holes) of 71, 72 independent functions (the partition function outside the region con- 
sidered above) times products of functions of the type ( |A.1| ) or ( |A.2| ), is symmetric in 71, 

72- 
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